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Abstract . This paper describes a flight-test methodology for developing a data base 
to be used to identify an aerodynamic model of a vertical and short takeoff and land- 
ing (V/STOL) fighter aircraft. The aircraft serves as a test bed at Ames for ongoing 
research in advanced V/STOL control and display concepts. The flight envelope to be 
modeled includes hover, transition to conventional flight and back to hover, STOL 
operation, and normal cruise. Although the aerodynamic model is highly nonlinear, it 
has been formulated to be linear in the parameters to be Identified. Motivation for 
the flight-test methodology advocated in this paper is based upon the choice of a 
linear least-squares method for model identification. The paper covers elements of 
the methodology from maneuver design to the completed data base. Major emphasis is 
placed on the use of state estimation with tracking data to ensure consistency among 
maneuver variables prior to their entry into the data base. The design and processing 
of a typical maneuver will be Illustrated. 


Keywords . Flight testing; data consistency analysis; aerodynamic modeling. 


INTRODUCTION 

Ames is conducting a flight research program on 
guidance, control, and display concepts for V/STOL 
aircraft. The program is directed towards extend- 
ing V/STOL capability to include flight operations 
aboard small ships in adverse weather. The test 
bed for this research is a YAV-8B aircraft, a 
prototype of the currently operational AV-8B Har- 
rier. The Harrier is a subsonic, vectored-thrust, 
V/STOL fighter aircraft; its engine nozzles can be 
rotated from 0° for forward flight to somewhat 
greater than 90° for hover and vertical flight. 
Compressor air is piped to the extremities of the 
aircraft for attitude control in hover and 
low-speed flight. The V/STOL Research Aircraft 
( VSRA ) is shown in Fig. 1. The purpose of this 
paper is to describe a flight-test methodology for 
developing a data base that will be used to iden- 
tify a full-envelope VSRA aerodynamic model. The 
model will be used to update and improve an exist- 
ing VSRA simulation {Anderson and Bunnell, 1985) 
which will aid in design of guidance, control, and 
display systems for the aircraft. 

The goal of acquiring a data base matched to a 
least-squares (regression) identification task has 
motivated the flight-test methodology outlined in 
Fig. 2. Because regression methods are computa- 
tionally simple, careful attention can be given to 
the structuring of an accurate and physically 
meaningful model. The VSRA data base consists of 
a set of individual 3-5 min flight-test maneuvers 
that together cover the complete flight envelope. 
The records from each maneuver are sufficient to 
construct time histories of the forces and moments 
that comprise the aerodynamic model, and include 
time histories of all relevant aircraft state and 


control variables. Long (15-30 min) records, 
consisting of concatenated segments from several 
maneuvers, will be used to identify parameters in 
each model equation. The important aspects of the 
preflight planning, flight testing, and postflight 
processing phases necessary to acquire the data 
base will be covered in this paper. The modeling 
of VSRA aerodynamics will be reported elsewhere. 

A key element in this methodology is the use of 
state estimation in the processing of each maneu- 
ver before its entry to the data base. Flight 
data often contain significant errors which should 
be corrected before the data are used for model 
identification. Furthermore, direct measurements 
of some dynamic variables may be unreliable or 
unavailable. These problems can be solved by the 
analytical method known as state estimation. 
Application of the method to aircraft flight data 
is possible because the forces and resulting 
motions of an aircraft along its flightpath are 
related by well-known equations of motion. The 
equations are used to produce estimates of force 
and motion variables that are compared with corre- 
sponding measurement time histories, usually with 
iterations, until suitable "fits" are obtained. 

The first application of state estimation to post- 
flight data analysis can probably be attributed to 
Otto Gerlach in the 1960s at the Delft Technologi- 
cal University, the Netherlands. This early con- 
tribution (Gerlach, 1966), called "flightpath 
reconstruction," was primarily concerned with 
accurate determination of angle of attack, pitch 
angle, and vehicle velocity during dynamic maneu- 
vers. These "states" were obtained by integrat- 
ing functions of measurements from the pitch-rate 
gyro and normal and longitudinal accelerometers. 
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Initial conditions and bias terms were determined 
from airspeed and altitude at the steady-state end 
points of the maneuver. The resulting "smoothed" 
time histories were then used as a basis for sub- 
sequent parameter- identification studies. 

In this country, Wingrove (Wlngrove, 1973) at NASA 
was an early advocate of state estimation for 
flightpath reconstruction. Over the past few 
years, work in this field has been evolving toward 
the use of more complete models, the development 
of more sophisticated algorithms, and the treat- 
ment of more difficult applications (Bach and 
Wingrove, 1985). State estimation methods are now 
used by many flight- test groups (Stalford and 
Ramachandran , 1978; Klein and Schiess, 1977; 
Hansen, 1979; Sri-Jayantha and Stengel, 1983). 

Once a consistent, smoothed set of time histories 
is obtained from the data, other analyses, such as 
aerodynamic model identification, are readily 
performed. The use of state estimation prior to 
modeling to correct the VSRA data and to provide 
estimates of unmeasured variables is illustrated 
in this paper. 

The paper proceeds as follows: the next section 

describes the planning and design of test maneu- 
vers needed to cover the flight envelope. The 
following sections discuss the data acquisition 
procedure used to generate the raw data records, 
and the techniques used to process those records. 
Then a section is presented that shows how state 
estimation is applied to the analysis of VSRA 
flight data. (The basic algorithm is given in the 
Appendix.) The last sections cover the calcula- 
tion of aerodynamic forces and moments, and illus- 
trate the procedure with a typical VSRA maneuver. 

PREFLIGHT PLANNING 

The VSRA Aerodynamic model must represent the 
three body forces and three moments over a flight 
envelope that includes hover, transition to for- 
ward flight and back to hover, as well as short 
takeoff and landing (STOL) operation and normal 
cruise. The resulting model will be strongly 
nonlinear with respect to aircraft state and con- 
trol variables such as angles of attack and side- 
slip, Mach number, nozzle angle, and power set- 
ting. The model nonlinearities can be conven- 
iently expressed with functions that are linear In 
the parameters to be identified (McNally, 1986; 
Anderson, 1986; Anderson and Hansen, 1987). A 
linear least-squares (regression) method (Klein 
and Batterson, 1983; Draper and Smith, 1981; Klein 
and Batterson, 1986) is well-suited to identify a 
highly nonlinear model that is linearly 
parameterized. 

The data base required for least-squares aerody- 
namic model identification can be obtained quite 
efficiently. Because the model is nonlinear, it 
is not necessary to maintain trim during a maneu- 
ver. In addition, because a regression procedure 
will be used to identify the model, large amounts 
of data may be processed. The order in which data 
samples are entered into the procedure is not sig- 
nificant, since an "equation-error" method does 
not require solution of the differential equations 
of aircraft motion. Accordingly, each flight-test 
maneuver has been designed to yield large changes 
in state and control variables while covering a 


(nearly) closed course within 5 min under continu- 
ous radar surveillance. The raw data base con- 
sists of as many longitudinal, lateral, and tran- 
sitional (to and from hover) maneuvers as are 
necessary to cover the flight envelope. After 
processing, model sections may be identified using 
long (15-30 min) records, each consisting of con- 
catenated segments from several maneuvers. 

One characteristic that sets the VSRA apart from 
conventional aircraft is that it exhibits signifi- 
cant power-induced aerodynamic effects when the 
nozzles are not in the full-aft position. These 
are largest during transition from hover to for- 
ward flight (and back to hover), and during 
periods of low-speed flight. Standard V/STOL 
procedures were used to provide data for identifi- 
cation of power-induced aerodynamics. One of 
these procedures, a short takeoff -slow landing 
maneuver, is outlined on the flight- test card 
shown in Fig. 3. In this maneuver, the ground 
roll begins with nozzles at 10°. At V r (indi- 
cated air speed) the nozzles are rotated to an 
angle e p (in the example for this paper, V r s 50 
knots and 8 p s 55°). Shortly after liftoff, the 
nozzles are rotated to the full-aft position. For 
the slow-landing portion, nozzles are rotated to 
40* Just before the final turn, and during the 
final approach are further rotated to 60°. 

Other maneuvers were specifically designed to 
excite large changes in longitudinal variables 
from several nominal trim points. In each of 
these maneuvers, nozzle angle, flap deflection, 
and power are held constant while the stabilator 
is varied to obtain changes in angle of attack and 
pitch rate. The maneuver includes "stick pump- 
ing," an "alpha ladder," and a "wind-up" turn. 

Near the end of the maneuver, power is added to 
return to the nominal trim point. Note that a 
considerable variation in Mach number may be 
experienced during the maneuver. Another set of 
maneuvers was designed to excite large changes in 
the lateral-model variables (angle of sideslip, 
yaw rate, and roll rate). All maneuvers were 
performed without "stability augmentation" to 
ensure a full range of vehicle response activ- 
ity. V/STOL procedures were performed in and out 
of ground effect. 


DATA ACQUISITION 

The VSRA measurement system is equipped with a 
10-bit digital data acquisition and telemetry (TM) 
system. A pulse-code modulation format is used to 
encode 156 mainframe channels sampled at 120 Hz 
and 160 subframe ‘channels sampled at 30 Hz. 

Before encoding, each analog channel is passed 
through a third-order Butterworth anti-aliasing 
filter with its cutoff frequency set at one-fifth 
of the channel sampling rate. After encoding, all 
flight data are transmitted to a ground station 
where they are recorded. A partial list of 
onboard measurements, those necessary for aerody- 
namic model identification, is given in Table 1. 

Flight tests of the VSRA were performed at the 
NASA test facility located at Crows Landing, 
California. The facility control room, which has 
a clear view of the runway and hover pad, is 
equipped with five eight-channel strip-chart 
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recorders and three color monitors for real-time 
display of the TH data. Two on-site radar systems 
are available to provide continuous tracking of 
test aircraft position. (A laser tracking system 
is used for all hover maneuvers.) During flight 
test, TH data from the VSRA onboard system are 
downlinked and merged at the facility with range, 
bearing, and elevation data from the tracking 
systems, and then recorded. 


PRELIMINARY PROCESSING 

Following real-time acquisition of data during 
flight test, each recorded maneuver, with instru- 
ment calibrations added, is converted to engineer- 
ing units and made available to researchers as a 
raw flight-data file. The first program in the 
post flight processing system is called PRODAT 
(PROcess DATa). This program reads the raw file 
and creates a "processed" file of selected chan- 
nels. Processing begins by removing wild points 
from the records. Several options are available, 
but one effective (and time-consuming) method is 
to pass each record through a "moving window." 
Points that fall outside the window are considered 
wild, and are tagged but not removed. When all 
wild points in a record have been tagged, the 
record is passed through a zero-phase-shift , low- 
pass digital filter to obtain an interpolated time 
history free of wild points. After interpolation, 
the data rate can be reduced to a submultiple of 
the mainframe sampling frequency. The filter 
cutoff frequency is set at one-half the final data 
rate desired. That data rate was chosen to be 
20 Hz for all VSRA maneuvers. 

Each channel processed from a maneuver raw-data 
file is stored in a processed flight-data file set 
up for that maneuver. The analyst may use a pro- 
gram called DSPDAT (DiSPlay DATa) to interactively 
select processed data channels for plotting in 
either x-y or strip-chart format. An x-y 
cross plot, for example, might display Mach number 
plotted against angle of attack. Such plots offer 
a convenient way to evaluate how well the flight 
envelope has been covered during a maneuver. It 
is unlikely that a single maneuver will provide 
enough variation in aircraft state and control 
variables to identify all model terms: the ana- 

lyst may also use this program to create a "map" 
file, which will contain addresses of time seg- 
ments selected from several processed maneuvers. 
This file can later be used to concatenate the 
selected segments to create a long record suitable 
for model identification. 


STATE ESTIMATION 

The next step in the processing of each maneuver 
is to apply SMACK (SMoothing for AirCraft Kinemat- 
ics), a state-estimation program developed at Ames 
Research Center (Bach, 1977; Bach, 1982) to check 
data consistency and derive unmeasured variables. 
The SMACK algorithm uses a six-degree-of- freedom 
kinematic model to fit aircraft measurement time 
histories. The algorithm is described in the 
Appendix. A diagram that defines the aircraft 
state and measurement models used in SMACK is 
shown in Fig. 4. Note that the state model is 
linear; all nonlinearities are found in the 


measurement model. Choice of a set of state 
variables and forcing functions consisting of 
Euler angles ( 4 , 8, 4), earth positions (x, y, h), 
earth winds (w x , w y , w h ), and their time 
derivatives leads co a linear state model. Note 
that the block symbols to the left of the vertical 
dashed line represent integrators. 

All nonlinearities associated with aircraft kine- 
matics appear in the measurement model. For 
example, the aircraft angular velocities (p, q, r) 
are expressed in terms of the states as 

p s i - 4 sin 8 

q s e cos 4 ♦ 4 sin 4 cos 8 (1) 

r s -0 sin 4+4 cos 4 cos 8 

Equations (1) are computed in the block of Fig. 4 
labeled with a script R. Calculations of air- 
craft body-axis velocities (u, v, w) and specific 
forces (a x , ay, a 2 ) require the usual transforma- 
tions from earth to body coordinates (Gainer and 
Hoffman, 1972) which are performed in the blocks 
labeled with a script L. Radar variables R 
(slant range), B (bearing angle), and E (eleva- 
tion angle), which are given by 

R . [(x - x a ) 2 ♦ (y - y a ) 2 ♦ <h - h a > 2 ] 1/2 

B . tan''[<y - y a >/(x - x a >] (2) 

E = sln' 1 t(h - h a )/R] 

where (x a , y a , h a ) is the antenna location, are 
calculated in the block labeled with a script P. 
If the measurement set includes static and total 
pressures, and total temperature, along with the 
air-flow angles, the measurement model can provide 
estimates of aerodynamic variables P (pressure 
ratio), V (airspeed), a (angle of attack), and 6 
(sideslip angle), from the relations (Dunlap and 
Porter, 1971) 

P * (1 - V 2 /7RT)~ 3,5 ; V = (u 2 + v 2 ♦ w 2 ) 1/2 

-1 -1 (3) 
a a tan (w/u) ; B * sin (v/V) 

where R is the gas constant and T is total 
temperature. These are calculated in the block 
labeled with script A. When air data measurements 
are included In the SMACK estimation procedure, 
the winds along the flightpath are also estimated. 
It is important that independent static pressure 
and air-flow angle corrections be available prior 
to a data-consistency analysis (Franklin, 1987). 


FORCE AND MOMENT CALCULATIONS 

The aerodynamic forces and moments acting on the 
VSRA during flight are determined as the differ- 
ence between the total forces and moments, and 
engine forces and moments. Here the term "engine" 
includes the reaction control system as well as 
the main nozzles. The engine forces and moments 
are calculated (offline) by a program called ENCAL 
(ENgine CALculations) . ENCAL uses a nominal pro- 
pulsion model of the VSRA Pegasus engine 
(YF402-RR-404) (Anon., 1983). Fan dynamics are 
not included in this version, since fan speed is 
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measured in flight. It should be noted that the 
propulsion model provides only thrust forces and 
moments. Any power-induced aerodynamic effects 
are to be included in the VSRA aerodynamic model. 

Inputs to the ENCAL routine include all the air- 
data, reaction-control, engine, and weight mea- 
surements listed in Table 1. Outputs to the pro- 
cessed flight-data file are the 3 body-axis compo- 
nents of engine forces and moments. ENCAL also 
calculates aircraft weight and inertias, and the 
variation in center-of-gravity location. These 
variables are added to the processed-data file. 

It should be emphasized that the aerodynamic model 
to be identified from flight data can only be as 
accurate as the engine model. A fully instru- 
mented Pegasus engine has recently been installed 
on the VSRA, and the engine model will be vali- 
dated after the next set of flight tests. 


Total VSRA force and moment time histories are 
obtained from the SMACK-derived estimates of 
accelerations and angular rates, and from ENCAL- 
derived estimates of weight and inertias. The 
body-axis forces are given by 


F s ma 
x x 


F y s ra y 


F .z 


( 4 ) 


In a final step, the aerodynamic forces and 
moments are calculated as the difference of total 
and engine forces and moments as outlined in the 
previous section. These are the time- histories 
that must be adequately represented by the VSRA 
aerodynamic model. Results of the ENCAL calcula- 
tions for the maneuver are shown in Fig. 7, with 
the corresponding aerodynamic variables shown in 
Fig. 8. Notice the tradeoff between engine and 
aerodynamic vertical forces during the STOL por- 
tions of the maneuver. 


CONCLUDING REMARKS 

This paper has described a flight-test methodology 
designed to generate a data base suitable for 
identifying an aerodynamic model using a 
regression procedure. Test maneuvers were 
designed to provide large changes in state and 
control variables around nominal trim points. The 
technique of state estimation was used to produce 
a consistent data set and to provide estimates of 
unmeasured variables. The availability of an 
accurate engine model has made it possible to 
isolate engine effects and determine aerodynamic 
forces and moments. 


where m is vehicle mass. The moments are calcu- 
lated from 


T l '• *xx a i - r zx (a n * - <I yy 


'zz*"' 


T „ » I yy a ra * W* ' P 2 > ' <‘zz ' < 5) 

T n s J z 2 a n " *zx^ a i “ “ ^xx " *yy^ M 

where I xx , I , I Z2 , I zx are vehicle moments of 
inertia. 


The completed VSRA data base consists of the set 
of individual 3-5 min flight-test maneuvers that 
covers the required flight envelope. Long 
(15-30 min) records, consisting of concatenated 
time segments from several maneuvers, can be 
created and used to identify parameters in each 
model section. An example of a typical maneuver 
in the data base has been shown in the paper. 

Work on developing the VSRA aerodynamic model from 
the data base is now in progress. 


EXAMPLE MANEUVER 

The short takeoff -slow landing maneuver described 
earlier serves to illustrate the type of informa- 
tion that is stored in the VSRA data base. In 
this maneuver, which contains abrupt changes in 
nozzle and flap angles, the aircraft transitions 
to normal flight after a short takeoff, performs a 
"go-around," and then transitions back to a STOL 
configuration for a slow landing. The raw data 
file includes all the onboard inertial and air- 
data, and radar tracking measurements indicated in 
Table 1. The variations in nozzle angle, flap 
setting, power, and control-surface positions 
required to perform the maneuver are shown in 
Fig. 5. Note how these time histories correlate 
with the activity requested of the pilot on the 
flight-test card of Fig. 3. 

The aircraft dynamic response to the control 
inputs is analyzed (along with the radar track) by 
the state-estimation procedure, which determines 
integrator initial conditions, selected instrument 
bias and scale-factor errors, and forcing-function 
time histories that provide the "best fits" to the 
measurement records (R, B, E), (P, o, 8), (a x , ay, 
a z ) , (p, q, r), and U, 8, 4>). The body angular 
accelerations (a t , a^ a p ), true airspeed (V), and 
flightpath winds (w x , w y , w h ) are also estimated. 
Results of the SMACK analysis required for calcu- 
lating forces and moments are shown in Fig. 6. 

4 


ACKNOWLEDGMENT 

The authors would like to thank programmers Mine 
Hagen and Suzi Kovacevich, and students Vafa 
Kordestanl and Phil Smith for their contributions 
to the work reported in this paper. We are grate- 
ful to Too Schultz for his work on the Pegasus 
engine model. 


APPENDIX: A STATE -ESTIMATION ALGORITHM 

In the context of this paper, "state estimation" 
refers to a process that solves a state model 

x = f(x,w) , x ( t Q ) = x Q (A. 1) 

such that h(x) in the measurement model 

z s h(x) + v (A. 2) 

suitably matches the data record over a time 
interval (t Q , t f ), usually in a least-squared 
error or minimum- variance sense. In Eq. (A.1), 
x is the state vector and w is a vector that 
represents any unknown forcing functions. In 
Eq. (A.2), h(x) is the system output vector; 2 and 
v are vectors representing the measurements and 
corresponding (random) measurement errors. 

For aircraft problems, the state and measurement 
models together represent the dynamics of a rigid 
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oody. The models are used to generate time histo- 
ries which include onboard variables such as Euler 
angles, angular rates, and linear accelerations, 
as well as tracking variables such as slant range, 
bearing, and elevation. Any bias or scale-factor 
errors associated with the state or measurement 
models are appended to the state vector and 
treated as constant but unknown parameters. 


The task of postflight state estimation prior to 
modeling is considered as a fixed-interval smooth- 
ing problem. Its solution consists of determining 
the x Q and w(t) that minimize the squared-error 
performance measure 


( x - xjV-V - x n )/2 

o 0 o o o 



<[z - h(x)] t R* 1 [z - h(x)] ♦ wV'w)dt/2 

(A. 3) 


subject to the dynamic constraint of Eq. (A.1). 

In Eq. (A. 3), x q is an a priori estimate of 
x Q , and P Q , Q? and R are weighting matrices. 
Mote that the first term of Eq. (A. 3) serves as a 
"penalty" function and tends to bias the estimate 
of x Q toward its a priori value. 


The necessary conditions for minimizing the per- 
formance measure lead to a nonlinear, two-point 
boundary value problem which is solved using a 
successive-linearization approach. A good algo- 
rithm involves choosing x Q and w, solving 
Eqs. (A.1) and (A. 2) to obtain a "nominal" 
solut ion, and then dete rmining a chang e in the 
initial condition from 




-( p :’ ♦ 

O O 0 0 


X ) ♦ a J 
o o 


(A. 4) 


where M 0 and a Q are found by solving (backwards 
in time) the "information” filter 


H . -f*M - Hf x - hjR-’h, ♦ 


f*a * hV’[z - h(x) ] ♦ Mf (Qf^a ♦ w) 

(A. 5) 


M(t f ) . 0 , a(t f ) . 0 


M 0 = M(t 0 ) , a 0 * a(t 0 ) 


In Eq. (A.5), M is the "information" matrix and 
a is the "gradient" vector, and f x , f w , and h x 
are Jacobian matrices. 


Next, a change in the forcing function is deter- 
mined by solving (forward in time) a first-order 
expansion of the dynamic constraint 

fix = f fix + f fiw , fix( t ) = fix 

X W o o (a6) 

with fiw - -w — Qf^(a + Mfix) 

The procedure is iterated until the performance 
measure of Eq. (A. 3) is minimized. It can be 
shown (Bach, 1977) that this algorithm is equiva- 
lent to the "modified” Newton -Raphson algorithm if 
there are no unknown forcing functions. 
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TABLE 1 Variable List for Aerodynamic Model Data 
Base 


Channel 

Measured 

Estimated 

Euler Angles 

Onboard 

SMACK 

Angular Rates 

Onboard 

SMACK 

Angular Accelerations 


SMACK 

Linear Accelerations 

Onboard 

SMACK 

Positions wrt Earth 

Radar 

SMACK 

Velocities wrt Earth 


SMACK 

Air-flow Angles 

Onboard 

SMACK 

Static Pressure 

Onboard 


Total Pressure 

Onboard 

SMACK 

Total Temperature 

Onboard 


True Airspeed 


SMACK 

Flightpath Winds 


SMACK 

Flap Setting 

Onboard 


Aileron Deflections 

Onboard 


Stabilator Deflection 

Onboard 


Rudder Deflection 

Onboard 


Engine Nozzle Angle 

Onboard 


Engine Fan Speed 

Onboard 


Compressor Pressure 

Onboard 


Fuel and Water Weights 

Onboard 


RCS Roll-Valve Positions 

Onboard 


RCS Pitch-Valve Positions 

Onboard 


RCS Yaw-Valve Position 

Onboard 


Engine and RCS Body Forces 


ENCAL 

Engine and RCS Moments 


ENCAL 

Gross Weight and Inertias 


ENCAL 



Fig. 1. The NASA V/STOL Research Aircraft (VSRA). 
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Fig. 3. Plan view of STOL Maneuver. 
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Fig. 2. Flow diagram for VSRA flight-test method- 
ology . 



Fig. 4. State and measurement models for the 
SMACK state-estimation procedure. 
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Fig. 5. Control tine histories for ST0L maneuver. 

(a) Nozzle, flap, power, (b) Stab i la tor, 
aileron, rudder. 
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Fig. 6. SMACK estimates for ST0L maneuver. 

(a) linear accelerations, (b) Angular 
accelerations. 
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Fig. 7. ENCAL estimates for STOL maneuver. 

(a) Engine forces, (b) Engine moments. 
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Fig. 8. Aerodynamic variables for STOL maneu- 
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